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characterized by the energy of the peak of vF^ and the evolution of this quantity 
is examined relative to the emission intensity. In most cases it is found that 
this peak energy either rises with or slightly precedes major intensity increases 
and softens for the remainder of the pulse. Inter-pulse emission is generally 
harder early in the burst. For bursts with multiple intensity pulses, later spikes 
tend to be softer than earlier ones indicating that the energy of the peak of 
vFy is bounded by an envelope which decays with time. Evidence is found that 
bursts in which the bulk of the flux comes well after the event which triggers 
the instrument tend to show less peak energy variability and are not as hard as 
several bursts in which the emission occurs promptly after the trigger. Several 
recently proposed burst models are examined in light of these results and no 
qualitative conflicts with the observations presented here are found. 

Subject headings: gamma rays: bursts — methods: data analysis 
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1. Introduction 

Gamma-ray burst continuum spectra provide the most direct information about 
the emission processes involved in these violent events. Unfortunately, the continuum 
generating process is unknown, making physical interpretation of spectral observations 
difficult. A striking feature of gamma-ray bursts is the temporal variability of spectra both 
between and within bursts. Therefore, an empirical study of the dynamics of burst spectra 
may provide vital clues for resolving the mystery surrounding these puzzling events. 

This is the second in a series of reports describing spectral observations of gamma-ray 
bursts as seen by the Burst and Transient Source Experiment (BATSE) on the Compton 
Gamma-Ray Observatory (CGRO). The goal of this series is to discover global properties 
of burst continuum spectra which may shed light on the gamma-ray burst problem. 
Gamma-ray bursts remain one of the least understood phenomena in astrophysics despite 
years of intense study. Their heterogenous nature makes classification difficult and few 
uniform burst properties are known. Bursts have only been observed above ;^2 keV and no 
emission has been seen in quiescence (|Schaefer 1994|) . In the first paper of this series ([Band 



et al. 1993| , hereafter Paper I), spectral properties of time integrated spectra from 54 bright 
bursts were studied. While spectra are quite diverse, it was shown that a simple and fiexible 
empirical model successfully described all spectra. This model is used in the present study 
to describe burst spectra and monitor the evolution of the energy at which the energy fiux 
per logarithmic energy band peaks [Ep, energy of the peak of uF^). In future work, similar 
evolution studies will be performed with other spectral characteristics, such as continuum 
shape or spectral bandwidth parameters. The data used here are high energy resolution 
spectra from the Spectroscopy Detectors (SDs) with moderate time resolution (> 0.128s). 
Subsequent work will use medium energy resolution data with fine time resolution from the 
BATSE Large Area Detectors (LADs) to examine weaker bursts and explore techniques 
using time-tagged counts from the SDs to improve time resolution. 

The study of burst spectral evolution has a long history. Golenetskii et al. (1983, 
hereafter G83) examined two-channel data covering ~ 40-700 keV with ~0.5 s time 
resolution from five bursts observed by the Konus experiment on Venera 13 and 14. 
Spectra were described by an optically thin thermal bremsstrahlung (OTTB) model, 
Ne{E) oc E~^ exp{—E / kT) ph-keV~^-s~^-cm~^. A correlation between luminosity and 
temperature parameter T was discovered (L ex T"', 7 ?«1.5-1.7), implying spectral hardness 
at a given time is related to the intensity of the burst in a simple way. Laros et al. (1985) 
performed a similar analysis using five bursts observed by the Pioneer Venus Orbiter 
(~ 0.1-2 MeV, ~0.2 s time resolution) but found no correlation and speculated that the 
G83 results were an artifact of the way temperature was inferred. 
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Norris et al. (1986, hereafter N86) investigated ten bursts seen by instruments on the 
Solar Maximum Mission satelhte using hardness ratios (the ratio of observed flux in two 
energy bands). The energy bands used in the hardness ratio were 52-182 keV and 300-350 
keV and the time resolution ranged over 0.128-ls, depending on the data available. They 
found that individual intensity pulses evolved from hard-to-soft with the hardness peaking 
before intensity. This implies a more complex relationship between hardness and intensity 
than proposed by G83. 

Recent studies of spectral evolution have tended to follow the lead of G83 or N86, 
obtaining similar results. Kargatis et al. (1994, hereafter K94) fit OTTB and thermal 
synchrotron models {Ne{E) oc exp[—(E/EcY^^] ph-keV~^-s^"'^-cm^^) to sixteen bursts 
from the SIGNE experiment, covering ~50-700 keV with 0.5 second resolution. A 
luminosity-temperature corelation was found in seven of the bursts (7 ~ 1.4-3.0) but 
two clearly did not show this correlation (the remainder were questionable) indicating 
that while many bursts may have a luminosity-temperature correlation, the trend is not 
universal. On the other hand. Band et al. (1992a) analyzed nine bursts observed by the 
BATSE SDs (~25-1000 keV, >0.128s) using models of the form Ne{E) oc E"^ exp[-{E/kT)] 
ph-keV~^-s~^-cm~^, confirming the results of N86 including T leading the intensity when 
fine time resolution was available. Bhat et al. (1994) used hardness ratios with BATSE 
LAD data to study single pulse bursts which had a fast rise followed by a gradual decay. 
The energy bands used for the ratio were ii^=25-100 keV and E >100 keV. The time 
resolution of the sample was 64 ms. Hard-to-soft spectral evolution was found with hardness 
leading the intensity. The time lag between the peak hardness ratio and peak count rate 
was found to be directly correlated with the rise time of the counting rate. 

Although the analyses exemplified by G83 and N86 give apparently inconsistent results, 
they can be reconciled by considering the time resolution of the data. The analysis of 
G83 and K94 had a time resolution of ~ 0.5 s while the resolution of N86 was as small as 
~ 0.128 s and Bhat et al. (1994) had 64 ms resolution. Therefore, slightly asynchronous 
behavior could be masked by the poor time resolution in the G83 and K94 samples. This 
conclusion is supported by the results of Band et al. (1992a), who performed an analysis 
similar to that of G83 and K94, but found that the spectral hardness (parameterized by T) 
led the intensity on short timescales. The need for fine time resolution is emphasized by 
Kouveliotou et al. (1992), who examined 22 bursts observed with the LADs and performed 
Fourier transforms on intensity profiles in different energy bands. The hard band was found 
to lead the soft by ~0.1 second. In addition, Kouveliotou et al. (1994a) reported significant 
variations in the hardness ratio at the 2 ms level in an extremely intense burst observed 
with BATSE. 
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The spectral analyses of G83 and N86 characterize the spectrum in very different ways. 
Both G83 and K94 use spectral models which are assumed to describe the gamma-ray 
burst continuum and whose parameters are interpreted as physically meaningful. Most 
often, the model chosen is one which cuts off rapidly at high energies (e.g., OTTB). It 
was demonstrated in Paper I that such models do not reproduce spectra observed by the 
BATSE SDs and conflict with the observation of high energy emission by other experiments 
( PVlatz et al. 1985| ; pchneid et al. 1992| ; [Hanlon et al. 19941) . Consequently, these models 



cannot be physically correct and though the parameters are certainly physically linked to 
the emission process, one must realize that the parameters do not necessarily have their 
apparent meaning (e.g., the OTTB cutoff energy may not be temperature). On the other 
hand, hardness ratios (such as used by N86) do not assume any knowledge of the source 
physics but do depend on instrumental properties which are accounted for in spectral 
fitting. While this characterization tracks changes in the shape and slope of a spectrum, the 
numerical values do not have a quantitative physical meaning and do not lend themselves 
to additional insights as a characteristic temperature or spectral shape can. Therefore it 
is helpful to use an empirical model which describes the photon spectrum but does not 
depend strongly on preconceived notions of spectrum formation in bursts. 

In this work, a description of the photon spectrum independent of physical models 
is achieved by fitting the empirical spectral model from Paper I to BATSE SD data. 
This model was shown to adequately describe spectra in the BATSE energy range and 
incorporates many simple physical models as special cases, accommodating uncertainty 
about the actual continuum emission process. Although the model is phenomenological, 
physically meaningful parameters can be derived from it. In this work, Ep (the energy of 
the peak of vFi,) is used to quantify fitted spectra, a parameter which indicates the energy 
of maximum radiated power. This allows a study of spectral evolution in gamma-ray 
bursts based on a well-defined physical measure of spectral hardness which does not depend 
strongly on assumed emission processes. 

In this work, the word hardness is used to describe the energy of Ep. Spectra in which 
Ep is at high energies are called hard and soft spectra have small Ep. Before discussing 
the analysis, the SDs are briefly described (§2). In §3 the analysis techniques used are 
presented in detail along with several consistency tests for both the data and methods. 
The results are presented in §4 followed by a discussion of their implications for both burst 
phenomenology and recently proposed theoretical models (§5). The entire work is then 
briefly summarized (§6). 
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2. The Instrument 

BATSE is a set of eight detector modules located on the corners of the CGRO 
spacecraft. Each module contains two detectors: an LAD, optimized for the detection and 
location of gamma-ray transients, and an SD, optimized for energy resolution. The SDs 
were designed for spectral studies and are used in this work because of their superior energy 
resolution to the LADs, which have large area but are thin (20" diameter by 0.5" thick) 
Nal(Tl) crystals. 

The gamma-ray detector for an SD is a 5" diameter by 3" thick Nal(Tl) crystal. 
Because of their thickness, the response to incident gamma-rays is roughly constant over 
a large range of energies and viewing angles. The SDs independently cover two energy 
decades in the range ~10 keV-100 MeV (the exact energy range is determined by a 
commandable gain setting, with higher gains covering lower energy ranges). The fractional 
full width at half maximum energy resolution of the SDs at 662 keV is ~7%, with a ~ E'^''^ 
dependence. A 3" beryllium window in the front of the aluminum case containing the Nal 
allows the response to extend down to ~5 keV for face-on bursts. In general however, the 
lower level discriminator is set so that the spectrum is cut-off below ~10 keV in the highest 
gain setting. An electronic artifact near the low end of spectra (the Spectroscopy detector 
Low Energy Distortion [SLED], Band et al. 1992b) exists in several channels above the 
instrument's low energy cutoff which can cause problems for spectral analysis. In this work, 
only data which are uncontaminated by this artifact are used. The energy deposited in the 
crystal is analyzed into 2782 linear channels which are rebinned on the spacecraft into 256 
quasi-logarithmic channels for transmission to Earth. Despite this compression, the width 
of most of the compressed channels is significantly smaller than the detector resolution and 
little information is lost in this process. 

BATSE enters a 4-10 minute long burst mode when the count rate in at least two LADs 
exceeds the background count rate by a predetermined significance threshold (usually 5.5(7) 
on one of three timescales: 64 ms, 256 ms, and 1024 ms. Several data types are collected in 
this mode. The type used in this study consists of high energy resolution SD spectra with 
accumulation time based on a time-to-spill criterion in the LADs (Spectroscopy detector. 
High Energy Resolution, Burst, called SHERB). When the total number of counts in the 
LADs exceeds a pre-determined threshold, the current SHERB accumulation ends and the 
next begins. A total of 192 spectra are recorded from the SDs associated with the four 
most brightly illuminated LADs. Half of these spectra are from the SD associated with the 
brightest LAD, one-quarter with the second brightest, and one-eighth each for the other 
two detectors. The length of the accumulation is a multiple of 64 ms with a minimum 
accumulation of 128 ms. Spectra are accumulated until either all spectra are exhausted 
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or the burst mode ends. The total time interval covered by each detector is equal: in the 
time required to accumulate one spectrum in the third and fourth brightest detector, four 
spectra from the brightest and two spectra from the second brightest detectors are recorded. 
Additional details about the BATSE detectors can be found in Fishman et al. (1989) and 
Horack (1991). 



3. Analysis 

It was shown in Paper I that the following empirical model is sufficient to describe 
burst spectra: 
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where A, a, (3, and Eq are fit to observed spectra and A' is chosen so that the function 
is continuously differentiable everywhere. The success of this model can be traced to its 
fiexibility. Included in equation (1) as special cases are a power law, photon exponential, 
OTTB, and a broken power law with continuous transition. Since the emission process is 
unknown, equation (1) is a logical choice for determining physically meaningful parameters. 
The quantity of interest in this study is the energy of the peak of vF^j (Ep) and since 
Ne{E)E^ oc uF^, Ep= (2 + a)Eo (assuming f3 < -2). 

Burst data were selected from detectors in high gain states, most covering ~ 15- 1500 keV 
with about a third extending to higher energies. This choice was motivated by the results of 
Paper I which demonstrated that for most bursts the spectral slope changes at energies less 
than 1 MeV and that the observed signal is greatest around 100 keV (a detector dependent 
property, not necessarily true for the incident spectrum). 

Bursts which had many high quality spectra were required so that Ep could be reliably 
determined and meaningful trends discovered. However, most bursts were not sufficiently 
long or intense to provide enough spectra to resolve the evolution. To understand the 
necessary signal strength, a Monte Carlo simulation of a typical spectrum was performed 
to investigate the reliability of the determination of Ep as a function of the signal-to-noise 
ratio (S/N) in the 60-200 keV energy band. The results of this simulation are given in 
Table 1 which shows that a large S/N level is required to reduce the variance in Ep to 
reasonable values. 
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Even for many strong bursts, most individual SHERB spectra do not have large S/N 
levels. However, the quality of the signal can be improved by averaging spectra together 
weighted by livetime. An explicit demonstration of this process for burst 2B910807 is given 
in Figure 1. In this figure, Ep as determined from spectra at the finest time resolution 
available and from averaged pairs of spectra are plotted together. S/N in the 60-200 keV 
band ranged from 7.0-17.0 with a median of 11.7 for the fine time resolution data and 
9.7-23.1 with median 16.3 for the averaged spectra. It is apparent that not only do the 
estimated errors in the determination of Ep get smaller as spectra are accumulated, but 
that the fitted value converges to the average of Ep determined from the fine resolution 
spectra. 

After selecting candidate bursts, background models were created using a channel-by- 
channel fit of a time dependent polynomial to data before and after the burst. The order 
of the polynomial varied from burst to burst; typically the lowest order which gave an 
acceptable fit was used (never higher than fourth order). After subtracting the background, 
consecutive spectra were averaged together until S/N>15 in the 60-200 keV energy band. 
Although Table 1 shows that this choice will result in a sizable dispersion in Ep, a larger 
S/N level would have essentially destroyed the time resolution in most of the sample. 
Therefore, some certainty in Ep was sacrificed for time resolution. Exceptions to this S/N 
criterion were made for a few bursts in which a long period of low activity was bounded 
by interesting pulse structure. In those cases, all data in between the intense periods were 
averaged despite a low S/N value. An exception was also made for burst 2B910813 in which 
the data did not extend below 100 keV. The S/N criterion for this burst was S/N>12 in the 
100-200 keV band. 

A burst was included in the sample if it survived the process outlined above with at 
least six spectra. Since most bursts do not have many spectra with large S/N, only ~4% of 
all bursts observed by BATSE and 17 of the 54 bursts from Paper I were contained in the 
sample (this sample includes bursts which occurred after those in Paper I). The empirical 
model (eqn. [1]) was fitted to data using forward folding deconvolution (P^oredo fc Epstein 



198^. This method determines the incident photon spectrum by folding an assumed 
photon spectrum through a model of the detector response and comparing the result to 
observed data. The best fit was found using a modified version of the Levenberg-Marquardt 
iterative x^ minimization algorithm ( [Bevington 1969| , p. 237; Press et al. 1989| , p. 521). 



The modifications were the use of model variances instead of data variances in x^ (see 
the appendix for details) and stopping criteria based on x^, x^ P^r degree of freedom in 
the fit {x^)-, and a test for false minima. The detector response model included the direct 
component along with scatter off the spacecraft ( [Pendleton et al. 1989D and the Earth's 



atmosphere ( [Pendleton et al. 1992|) . The error in Ep was determined through standard 



error propagation using the error in a and Eq and the correlation between these two 
parameters ( [Bevington 1969| , p. 61). 



Each spectrum was fit over the entire energy range above the SLED. The data were 
not rebinned from the telemetered format so that in most cases the average number of 
counts per bin in the high energy channels was less than a few. In these channels, the 
Poisson statistical distribution is not well approximated by a normal distribution and 
the x^ statistic is therefore not strictly applicable. This is a common difficulty in high 
energy astrophysics which can be avoided by binning the data until the Gaussian limit is 
reached. To determine the importance of this problem, a set of Monte Carlo simulations 
was performed for several types of binning. The results for four different binning modes in 
Table 1 show that the fits are not improved by rebinning the data even if most high energy 
bins have only a few counts. There are two reasons for this. First, the data at higher 
energies are so sparse that when the Gaussian limit is achieved, the bins are so wide that 
the upper power law in eqn. (1) cannot be accurately determined. Second, the fit is largely 
determined by the signal below a few hundred keV where the number of counts per channel 
is largest (Max(C) in Table 1). The Gaussian limit can also be achieved if there are a large 
number of background counts, which can be important at low energies. 

These simulations demonstrate that for the SHERB data type used in this sample, any 
binning is valid. To ensure that this was the case for real data, fits using the telemetered 
binning and eight broad energy bins were made to spectra from burst 3B930916. The 
comparison of the values of Ep is presented in Figure 2. Ep is consistent between the 
two fits and determined with roughly equal accuracy. Since the choice of data binning is 
unimportant, the telemetered channels were not rebinned because the analysis software was 
optimized for this scheme. 

One difficulty with using high gain detectors (where the spectrum extends to low 
energy) was that occasionally the upper spectral index /? was larger than —2 which implies 
that Ep was beyond the upper end of the fit range, that the spectrum has more gradual 
curvature than the empirical model could accomodate, or that the high energy signal was 
insufficient to fix /5. For this work, /3 was constrained to be less than —2 for all fits. In 
bursts for which this constraint was important (implying a fiatter spectrum), it was found 
that either /3 tended to shift rapidly from the upper constraint to the lower (/3 = —5) or 
that it remained at the upper constraint (—2.01) consistently. In the former case, the high 
energy spectrum was probably not well-determined and should be ignored when calculating 
Ep. For the latter, Ep should refiect changes in the shape (i.e. hardness) of the spectrum 
since Eq will move to higher energies to compensate for the inability of /5 to adequately 
describe this region. No bursts were observed to show convincing evidence for evolution of 
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13 from above —2 to lower values. Therefore, variations in the fitted value of Ep should be 
consistent with the true variations in spectral hardness if quantitatively inaccurate (with 
systematically smaller values of Ep). Omitting the upper power law would force Ep to 
higher energies as the model accounts for power at the high end. However, this simpler 
model may overestimate Ep, placing it far above the detector's energy range where the 
model is not as sensitive to changes in Ep. As a subset of the empirical four parameter 
model of equation (1), the model without j3 is not a more accurate description of the 
true photon spectrum. Since evolution of Ep is the primary focus of this study, equation 
(1) with (3 constrained is used to model all spectra since it remains sensitive to hardness 
variations over a broader energy range than a model which cuts off at high energies. 

To test whether changes in the fitted values Ep were indeed consistent with real 
variations in Ep for cases where 13 > —2, detectors at lower gain were used to determine 
Ep for bursts in which the constraint was important. Unfortunately, in all but one case 
the high energy signal was too weak to reliably fit Ep with any time resolution. Figure 3 
shows a comparison of the high and low gain detectors for the exception (2B9 10503). It is 
apparent from this figure that not only do the two detectors have similar time histories, but 
that Ep determined by the high gain detector is not far from the value of Ep measured 
by the low gain detector (for which (3 < —2). Therefore, it is not expected that qualitative 
statements about the evolution of Ep depend strongly on the accuracy with which (3 is 
determined. 

That the low gain detector measures a consistently larger value of Ep than the high 
gain detector in the previous example suggests the calculated value of Ep may depend on 
gain setting even though Ep is well defined in all detectors. To ensure that this was not 
the case, the value of Ep from two detectors at different gains viewing burst 2B921207 
were compared {(3 < —2 for both detectors). In this case, spectra accumulated by the high 
gain detector were averaged so that the time intervals covered by both sets of spectra were 
identical. The comparison in Figure 4 shows that the determination of Ep was consistent. 

The tests outlined above cover the obvious difficulties involved in the analysis. 
Although solutions to the S/N and f3 > —2 problems are not ideal, they do allow the 
analysis to proceed so that useful information about burst continua can be derived. Also, 
the inter-detector test indicates that SD observations are consistent with each other and 
the results should not depend on the variable energy range covered by the SDs. 



Results 
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Of the 862 gamma-ray bursts detected by BATSE until the end of 1993, only the 37 
bursts listed in Table 2 met the criteria outlined above. As expected, equation (1) fit the 
data well, xl ^ ^ for nearly all spectra. For these bursts, Ep was compared to the count 
rate (see Figure 5). Two bursts which clearly demonstrate the general trends observed 
in the sample as a whole are 2B921207 and 2B920525. In 2B921207, it is apparent that 
Ep softens over the whole burst except for the increase along with the intensity 8 seconds 
after the trigger. 2B920525 is a burst in which two strong intensity pulses were resolved. 
Ep softens after intial hardening within both spikes and the second spike is softer than the 
first. 

The plots of Figure 5 were used to characterize types of evolution in Ep. The 
classification categories used were motivated by the work of G83 and N86 in which a 
spectral characteristic (i.e., OTTB temperature and hardness ratio respectively) was tied 
to changes in burst intensity. These categories describe both the tendency for variations in 
Ep to refiect changes in intensity (G83) and hard-to-soft evolution (N86). The categories 
are: 

1. An increase in Ep occurring in proximity to a major intensity increase. Ep and 
intensity need not be morphologically identical and Ep can lead or lag the intensity 
increase by a small amount (^ 1 s). 

2. General softening of Ep in time outside of intensity pulses over the entire burst. 
Regions which do not have a spike-like intensity profile are compared for this category. 
Spikes are considered to be short time intervals (less than a few seconds) in which the 
count rate was much greater than the rate in nearby intervals. For broader pulses, the 
first few seconds after the maximum count rate was achieved are considered part of 
the pulse, the remainder was considered inter-pulse emission. 

3. Softening of Ep within intensity pulse structures. Ep should be harder early in the 
intensity pulse for a burst to fall in this category. It can be measured only if the pulse 
is temporally resolved. 

4. Later intensity pulses have a softer peak in Ep than earlier pulses. This can only be 
measured if the burst has multiple intensity spikes. 

Bursts were judged as showing the property in question, showing an opposite trend, or 
showing neither trend (+, — , and respectively in Table 2). To demonstrate how a 
burst was categorized, consider 2B921207. This burst satisfies the first category because 
Ep hardens as the intensity increases at 0.5 and 8 seconds after the trigger. The interpulse 
region lies between 2.5 and 8 seconds, where Ep clearly softens. Both intensity spikes 
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(0.5-2.5 seconds and 8-10 seconds after trigger) are temporally resolved and Ep softens in 
each. Finally, the peak Ep is softer in the later spike than in the first. 

The rating system used here is admittedly simplistic, but is preferable to schemes 
which allow for finer distinctions within a class (e.g., rating bursts on a scale of 1-10, 
either discretely or continuously) for two reasons. First, the sample size is small so that 
the subdivision of bursts into several discrete subcategories reduces the significance of the 
classification. Second, the mediocre time resolution for some time intervals and sizable 
errors for some Ep create problems for a continuous classification scheme. Therefore, 
although the criteria outlined above are subjective, they are the least problematic given the 
state of the observations. 

Individual bursts are detailed in Table 2 and characteristics for the sample as a whole 
are given in Table 3. Included in Table 2 are the median values of Ep and the range in 
Ep [AEp=m.ax{Ep)-m.m{Ep)] for each burst. Not all bursts could be classified because of 
large errors in Ep, unresolved structure, lack of structure, poor time resolution, etc. Bursts 
in which (3 was fit consistently at the upper constraint are noted in Table 2 as are bursts 
which were simultaneously observed in the 1-30 MeV range with the COMPTEL instrument 
on CGRO ( lianlon et al. 19"U^ ). In all four cases, the fitted (3 was consistent with power law 



fits to the COMPTEL data in that for bursts in which COMPTEL reported a power law 
of index > —2, (3 was consistently fit to the upper constraint. It can be seen from Figure 5 
that there are several instances in which Ep was poorly fit resulting in a large error {(Tep) 
for this quantity. Spectra for which aEp ^0.5Ep were ignored when characterizing bursts 
and calculating AEp. 

Table 2 clearly demonstrates that increases in Ep are associated with rises in intensity 
and examination of Figure 5 shows that in many cases, the peak in Ep leads the intensity 
peak. A quantitative analysis of Ep leading intensity was not performed since Kouveliotou 
et al. (1992) found the typical lead time is shorter then the resolution of this sample. In 
addition, the time profile is extremely complex on short timescales so that a detailed study 
would best be performed using data with better time resolution. This data can be obtained 
from the LADs and the analysis is currently being performed (Kouveliotou et al. 1994b). 
Table 3 shows that many bursts undergo hard-to-soft evolution both within and outside 
of intensity spikes as well as between successive pulses. Softening of successive spikes is 
explicitly demonstrated in Figure 6, where the peak in Ep for early intensity spikes is 
compared to the same quantity in later spikes in multi-spike bursts. It can be seen that 
late pulses are more likely to be softer than those which occurred previously, suggesting 
that the maximum possible value for Ep is bounded by an envelope decaying in time. The 
trends noted above suggest that spectral hardness leads the intensity modulated by a decay 
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envelope. 

It should be mentioned that there are bursts in the sample which do not share the 
traits outlined above. This shows that while typical burst characteristics may have been 
identified, they are not universal on the timescale at which this sample was studied. 
However, it is not unreasonable to expect that Ep leading intensity and the softening trends 
noted above will be found in most bursts when observed at finer time resolution. 

The general softening trends in Ep imply that the time at which emission occurs 
relative to the beginning of the burst is an important factor in spectral evolution. This 
raises an interesting question: assuming that the maximum values of Ep are bounded by 
an envelope decaying in time, then if a burst has a precursor several seconds before most 
of the emission occurs, does this envelope begin to decay when the precursor happens or 
is the precursor irrelevant with respect to spectral evolution? To examine this question, 
the range in Ep and the time at which significant emission occurred relative to the BATSE 
trigger were compared (Figure 7). Significant emission was defined as emission capable of 
producing spectra which met the S/N criterion outlined in §3. In Figure 7 the horizontal 
bar indicates the period of significant emission. It begins at td, the time delay between the 
instrument trigger and the beginning of significant emission, and has length equal to the 
duration of significant emission (AT). A large value of td results when a long period of 
time passes between the trigger (a burst precursor) and the brightest emission. The data in 
Figure 7 occupy a triangular region indicating the sample was devoid of bursts with hard, 
highly evolving spectra which began emitting significantly well after the BATSE trigger: the 
few bursts with large td have small values of Ep and AEp. To determine the importance of 
these bursts, the sample was divided into four groups based on the values of td and AEp. 
Since there are no clear groupings in Figure 7, the divisions were t^ = (significant emission 
at trigger) and AEp=400 keV, which each divide the sample in half, producing an unbiased 
division. As shown by Table 4, bursts with t^^ 7^ do tend to evolve less, but the significance 
of this result is small (~ la), implying that the apparent td—AEp association may be a 
statistical aberration which will vanish when more bursts with large td are observed. 

For several bursts, the nonzero time delay was short compared to the duration of 
significant emission. This suggests that the measure of time delay should be the ratio 
of td and At. This definition is preferable to the absolute measure since it incorporates 
the timescale of the event from the trigger to the end of significant emission, allowing 
comparisons based on the time the burst began its most active stages and the total length of 
the burst. Bursts were divided into two groups based on the value of td/At. Table 5 shows 
the number of bursts in each group for several values of this ratio as well as the median 
and extreme values of AEp for each group. The table shows that late emitters tend to 



-13- 



show less Ep variability than bursts which emit promptly after the trigger. Although this 
is very suggestive, there is sizable overlap between the groups. Considering the small size of 
the sample, this result is not too significant and therefore the possibility of a dependence of 
AEp on td should be treated with some caution. 



5. Discussion 



In this work it was shown that for most gamma-ray bursts, the energy of the peak of 
uFi, is associated with and sometimes leads the intensity of the burst modulated by an 
envelope decaying in time. Indications that bursts whose main emission occurs at relatively 
long times after the trigger are generally soft and tend to evolve less in Ep were also found. 
In this section the phenomenological implications of these results are discussed. Some 
recently proposed burst models are also examined in light of the observations reported here 
and no qualitative conflicts between the models considered and observations are found. 



5.1. Phenomenological Implications 

While it has been known for some time that burst spectra tend to evolve from hard to 
soft (N86), the observations presented here show that the time of emission relative to the 
beginning of the burst may also determine how the spectral hardness evolves. For optically 
thin burst models (in which particles radiate on timescales much shorter than the resolution 
of SHERB spectra), these results imply that as the burst transpires, the emission region 
'remembers' previous emission suggesting that either emission comes from only one region 
or that different emission regions must be physically connected. The softening of successive 
intensity pulses could be caused by significant change in the emitting region during intense 
episodes (such as an expansion or an increase in the number density of particles). The 
softness of delayed emission implies that the emission region can evolve without observable 
radiation. For optically thick models, a burst could be the result of a single energy input 
since the radiation does not escape the region on short timescales. If the region cools by 
adiabatic expansion and the surface of last scattering does not propagate inward to higher 
temperature regions too quickly, then observable radiation would be expected to evolve 
from hard-to-soft. The softness of delayed emission could also be explained in this way if 
the region was not able to emit significantly during the initial stages of expansion. 
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The hard-to-soft spectral evolution within pulse structures indicates that intensity 
spikes are not symmetric in time. Nemiroff et al. (1994b) reached a similar conclusion by 
examining the time structure of intensity profiles. These observations clearly show that 
time symmetric models of spectra formation as well as overall intensity profiles cannot be 
correct. An example of the kind of model eliminated by these observations is symmetric 
beam models in which the beam sweeps by much faster than the timescale over which the 
beam changes. 

Another implication of this work relates to burst duration classes. The distribution of 
gamma-ray burst durations observed by BATSE is bimodal and short bursts (duration ^ 2 s) 
tend to have larger hardness ratios than longer ones ( [Kouveliotou et al. 1993| ; [Lamb et al.| 



199^ where the hardness ratios were averaged over the whole burst. Based on the large 
values of AEp reported here and the fact that the hardest emission tends to come early 
in a burst, short bursts are expected to be harder than long ones since they do not have 
softer tails, not necessarily because of an inherent difference in the radiation mechanism. 
An example of the effect of time averaging is burst 2B920525. If spectra from 4.16-6.72 s 
(the first spike) are used ii^p=650±70 keV but falls to 421±25 keV if the entire burst is 
integrated into one spectrum. As Norris et al. (1994) also suggested, if hardnesses are 
measured over the same absolute time interval (e.g., compare short bursts to the first few 
seconds of long bursts), the hardness distribution may be the same for both classes. 

In Paper I it was shown that Ep varied considerably from burst to burst indicating 
that there is no universal characteristic energy of burst emission. One conclusion which 
was drawn from this was that pair processes could be directly observed only if burst 
sources had a broad range of redshifts. Such processes would produce spectra with an 
annihilation feature at 511 keV and would cut off above this energy because of pair opacity. 
The observed energy of these features depends on the redshift of the source region. The 
argument against observable pair processes can be extended by noting the large changes in 
AEp observed within many bursts in this sample (Table 3). Unless the emission process 
involves large and rapid changes in redshift, pair processes cannot be directly observed in 
burst spectra. 

The data in Figure 7 include several bursts in which the emission was prompt but had 
small AEp. An interesting possibility is that these prompt emitting, relatively non-evolving 
bursts were preceded by a precursor too faint to trigger BATSE. In this case, bursts which 
occupy the lower left-hand corner of the figure are late emitters for which BATSE did 
not trigger on the initial event. A less exciting possibility is that bursts which emit late 
but evolve considerably exist and there is no relation between emission time and AEp. 
Although this type of behavior is not seen, it cannot be ruled out because of the small 
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sample size and rarity of late emitters (about one strong burst per year with td >60s). 
Therefore, this may remain an open question even if BATSE detects many bursts over the 
next few years. 



5.2. Implications for Burst Models 



A plethora of gamma-ray burst scenarios have been proposed ( Nemiroff 1994a| ) 



However, most attempt to explain only general spectral characteristics and energetics with 
little or no mention of how spectra are expected to evolve. Therefore, several models 
are examined to find those which can accommodate the observations presented here. 
Because BATSE observations of burst isotropy and inhomogeneity have cast serious doubt 
on older models, only those proposed after 1991 (the year BATSE was launched) are 
considered. Discussion is further restricted to models which purport to explain all or most 
gamma-ray bursts. Although arguments have been made for two population distributions 
( [Lingenfelter fc Higdon 19^ ; [Smith fc Lamb 1993) ) , the inability to separate bursts into two 



uncontroversial classes makes such discussion premature. Jet models are also avoided (c.f. 
Brainerd 1992| ; |Dermer fc Schlickeiser 1994| ) since the phenomenon is poorly understood 



and the models are incomplete. 'Exotic' sources such as cosmic strings ( Paczyhski 1988 



strange stars ( Paensel et al. 1991| ), or primordial black holes QCline fc Hong 1992|) are also 
omitted since the existence of such sources is highly speculative. 

The distribution of bursts seen by BATSE has forced bursts out to extended halo or 
cosmological distances if they arise from a single source population. Since source regions 
are inferred to be small from the fast rise times of bursts, the energy density at the source 
must be enormous, particularly if the burst radiates isotropically. Therefore, many models 
incorporate optically thick relativistic electron-positron fireballs. As originally proposed 
( poodman 1986| ; Paczyhski 1986|) , fireballs cannot explain the observed nonthermal emission 
and will not be relativistic if the fireball is contaminated by baryons ( Pavallo fc Rees 1978 ; 
Paczyhski 1990|) . 



Many different scenarios have been proposed recently (c.f. [Narayan, Paczyhski, fc 



Piran 1992| ; [Usov 1992t [Woosley 1992|) which can essentially be reduced to a fireball which 



must either form in a region free of baryons or break through a cloud of baryons to be 
observed. In these cases, the observed nonthermal spectra are produced either by several 
thermal blackbodies at different temperatures (corresponding to different fireballs) or 
through interactions with a strong or turbulent magnetic field. If a number of fireballs 
occur, then hard-to-soft spectral evolution could be explained by the creation of low 
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temperature fireballs or rapid cooling as the burst transpires. For magnetic interaction 
models, the general softening of Ep may reflect a gradual decline in available energy as the 
field evolves into a stable configuration. 

Meszaros & Rees (1992) have proposed a model involving a baryon contaminated 
fireball from any source, galactic or cosmological. The baryons in the fireball interact with 
the interstellar medium (ISM) to form a relativistic shock where gamma-rays are radiated. 
Hard-to-soft evolution could be explained either by the most energetic baryons reaching the 
shock first or the gradual decline of average energy as more particles encounter the shock. 
The softness of late emission could be a consequence of the baryon front losing energy to 
the ISM before the shock is fully formed or could arise from the shock losing energy as it 
propagates through the ISM. 

In a different type of model suggested by Melia & Fatuzzo (1992), a gamma-ray burst 
is generated by sheared Alfven waves from a radio pulsar. The waves are created by 
crustal disturbances at the polar cap which flood the magnetosphere with charged particles. 
These charges upscatter ambient radio photons to gamma-ray energies which are beamed 
outward. Hard-to-soft evolution may result from the most energetic charges scattering first 
and successive spikes soften if subsequent disturbances are less energetic than the original 
disturbance. Since the emission is beamed, the hard emission might occur when the beam 
was not pointed in our direction. Delayed emission might be soft because the early hard 
emission was not beamed in our direction. 



6. Summary 

In this work, the evolution of the energy of the maximum flux per logarithmic energy 
band (Ep) in long, bright gamma-ray bursts was studied using high energy resolution 
spectra from the BATSE SDs. It was found that Ep is coupled with the intensity of a burst 
(either leading or accompanying intensity increases) but is apparently modulated by an 
envelope decaying in time. These results are consistent with previous studies by N86 and 
Bhat et al. (1994). Indications were also found that bursts in which significant emission is 
substantially delayed relative to the instrument trigger tend to be soft and evolve less in 
Ep although the significance of this result is low. Several burst models were examined in 
light of these results and it was found that none of those discussed could be ruled out on 
the basis of these results alone. However, as the models mature to the point of more precise 
spectral predictions, the observations presented here may prove to be a useful constraint or 
suggest directions for further effort. 
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A. Model Weights for x^ 

The Levenberg-Marquardt algorithm (Bevington 1969, p. 237; Press et al. 1989, p. 
521) is a prescription to efficiently minimize goodness-of-fit statistics. The most common 
statistic used is 



x^ = E 



N 

Vi - y[Xi; a) 



(Ti 



(Al) 



where: x is the independent variable (here energy); y is an observed variable (count rate 
in this case); y{xi; a) is the value of the model at Xi for parameters a; A^ is the number of 
measurements (here, channels in a spectrum); a represents the parameters to be fit; and 
ai are the errors associated with yi. x^ is minimized using an iterative gradient procedure 
for which the first and second partial derivatives of x^ with respect to a are needed. In its 
usual implementation, the algorithm assumes that the a^ are taken directly from the data 
and are constant. 

The problem with using data variances (o"j) in equation (Al) is that ai weights 
downward fluctuations in the data too strongly, a result of the Poisson nature of count 
data in which af oc Ui, the number of observed counts. In this case, a small number of 
counts create a smaller than average o"j which causes the associated data points to be more 
important in the calculation of x^- Wheaton et al. (1994, hereafter W94) examined this 
problem in detail using a Monte Carlo simulation of 24,000 one-parameter data sets and 
found that the mean of fitted values deviated from the true value by —170(7 (note the 
downward shift). The solution to this problem is to use ai generated from the model under 
consideration rather than the data, replacing af ex rii with af{xi;a) oc ni{xi;a) (model 
variance). When W94 performed the simulation with this new scheme, the bias was — 0.34cr. 

Although W94 demonstrated the importance of using weights derived from the model, 
they do not mention a necessary change in the standard Levenberg-Marquardt prescription 
to accommodate the new weighting. The algorithm assumes that ctj is taken from the data 
so that changes in the x^ surface are contained in only the y{xi; a) term in equation (Al). 
However, by making o"j model dependent without changing the partial derivatives, assumed 
changes in the x^ surface caused by variation in a will be incorrect, resulting in slightly 
off-course iterations and an incorrect minimum since dx'^/da=0 can be at the wrong value 



-18- 



of a. It was found that when this effect is not accounted for, the true minimum of the x^ 
surface is not reached and the values of the fitted parameters differ from the true best fit 
by roughly the quoted error in the fitted parameters though it is not biased toward higher 
or lower values making it difficult to detect in large collections of fits. 

Alterations in the second derivative terms only affect the route used to find the 
minimum of x^ and therefore need not be corrected for the new weighting scheme ( [Press et" 



al. [1989| , p. 523). On the other hand, the first partial derivatives need to reflect any changes 



which might occur in the x^ surface and must be corrected. The corrected version is 



daj 



N 



-2E 



Vt 



af{xi;a) 



l' dyi{xi;a) _ [yj - yj^Xj; a)] daijxi; a) 
\ daj ai{xi;a) daj 



(A2) 



The first partial derivative term in this equation is equal to the partial derivative specified 
by the Levenberg-Marquardt algorithm with model variances replacing data variances but 
the second term is new. For this work, equation (A2) can be expressed in a simpler form 
since background subtracted count rates are used (y.j = rii/ fi — hi where /j is a term which 
converts counts to a count rate and hi is the background counting rate). In this case, 
crfi^Xi] a) becomes 

a^[Xi]a) = + 0-^ (A3) 

Ji 

where af. is the error in the determination of the background rate in channel i. 
Differentiating (A3) with respect to a and substituting into equation (A2), it is found that 

daj ^ cxf(xi;a) daj 

This expression is equal to the uncorrected expression when ctj (xj ; a) is replaced by o", 
which allows the standard algorithms (e.g. Bevington 1969, p. 237; Press et al. 1989, p. 
521) to be modified easily. When this expression is substituted into the algorithm, the true 
minimum of the x^ surface with model variances is found. 
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Table 1. Monte Carlo Simulations 



Binning^ 



tis) 



Med(C)3 Max(C)4 Med{Epf Yd.i{Epf S/N^ 



nobin 



0.1 
0.25 
0.5 
1.0 
2.0 



0.23 
0.99 
2.11 
4.41 
9.70 



4.66 

7.78 

16.48 

30.31 

43.72 



259 
285 
265 
251 
250 



286 
77 
50 
36 
21 



10.0 
15.6 
21.6 
29.6 
39.5 



half 



0.1 
0.25 
0.5 
1.0 
2.0 



1.18 
2.78 
4.66 
9.05 
21.50 



9.52 

26.85 

45.71 

86.85 

219.33 



271 
259 
256 
252 

251 



176 

85 

56 

39 

22 



10.0 
15.6 
21.6 
29.6 
39.5 



full 



0.1 
0.25 
0.5 
1.0 
2.0 



2.46 
4.95 
8.55 
19.30 

45.47 



15.47 

36.03 

72.05 

141.77 

321.04 



252 
244 
253 
251 

250 



153 

68 

50 

35 

21 



10.0 
15.6 
21.6 
29.6 
39.5 



eight 



0.1 
0.25 
0.5 
1.0 
2.0 



3.29 

8.54 
17.67 
34.85 
73.27 



10.52 

24.35 

46.76 

102.28 

195.23 



220 
240 
247 

250 

248 



86 
63 
46 
40 
25 



10.0 
15.6 
21.6 
29.6 
39.5 



^The binning scheme used for spectra: 

nobin — Data in native format. 

half — Binned to one half of the detector energy resolution. 

full — Binned to the detector energy resolution. 

eight — Eight broad energy bins. 

^The livetime used to create spectra. 

^Median number of counts per channel in the binned spectrum. 

■^Maximum number of counts per channel in the binned spectrum. 

^Median value of Ep (keV). 

^Variance of Ep about the median (keV). 

^Median S/N in the 60-200 keV energy band. 

Note. — 200 spectra were created using the model spectrum of eqn. (1) with 
A=0.1 ph-keV-^-s-i-cm-2, « = -1, /5 = -3, Eq = 250 keV {Ep=250 keV). Non- 
integral counts result from subtraction of fractional background counts. 
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Table 2. Classification of the Sample 



Burst 


1 


Characteristic 
2 3 


4 


Med{Ep) 
(keV) 


AEp 
(keV) 


2B910503*t 


+ 


+ 


+ 


+ 


540 


1350 


2B91060lt 


+ 









700 


950 


2B910807 


+ 





+ 





170 


175 


2B910814 










490 


900 


2B910814Ct 




+ 






1400 


1300 


2B91103r 




+ 






400 


850 


2B911106 











165 


120 


2B911118t 


+ 


+ 


+ 


+ 


180 


440 


2B911126* 


+ 


+ 


+ 





250 


900 


2B911127 










190 


250 


2B911202* 





+ 







350 


1300 


2B911209* 


+ 


+ 


+ 


+ 


270 


400 


2B920218* 




+ 






220 


150 


2B920311 


+ 


+ 


+ 


+ 


440 


1150 


2B920406* 


+ 





+ 





240 


300 


2B920513* 


+ 




+ 




220 


200 


2B920517* 


+ 




+ 




230 


150 


2B920525* 


+ 


+ 


+ 


+ 


380 


2075 


2B920617 




+ 






145 


175 


2B920622 


+ 





+ 




440 


1550 


2B920627* 










210 


325 


2B920711 


+ 










510 


650 


2B920720 


+ 


— 


— 


— 


270 


525 


2B920902* 




+ 


+ 




550 


1500 


2B921009 


+ 









300 


225 


2B921123 


+ 


+ 






250 


350 


2B921207 


+ 


+ 


+ 


+ 


140 


650 


2B921209 


+ 


+ 


+ 


+ 


180 


375 


2B930120 










125 


80 


3B930405 




— 






280 


300 


3B930425 










250 


225 


3B930506 


+ 


+ 







1030 


1750 


3B930916* 














390 


350 


3B930922 


+ 


+ 


+ 





110 


130 
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Table 2 — Continued 



Characteristic Med{Ep) AEp 

Burst 12 3 4 (keV) (keV) 



3B931103 


+ 


+ 


+ 


+ 


380 


1300 


3B931126 


+ 


+ 


+ 


+ 


240 


300 


3B931204 


+ 


+ 


+ 





380 


775 



*j3 = —2.01 throughout the burst. 

t Observed by COMPTEL ( |Hanlon et al. 1994|) 
Burst Characteristics: 

1 — ii^p-intensity association. 

2 — Overall softening within a burst. 

3 — Softening within pulse structures. 

4 — Later spikes softer than earlier spikes. 

The symbols denote: 

H Trait observed. 

Opposite behavior observed. 

— Neither trait nor opposite observed. 
... — Could not be determined. 
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Table 3. Summary of Classification 37 Bursts 



# Observed # Possible Property 

23 25 Ep- intensity association 

Ep softens over whole burst 
Ep hardens over whole burst 

Ep softens within intensity spikes 
Ep hardens within intensity spikes 

Later spikes softer than earlier ones 
Later spikes harder than earlier ones 



20 


30 


2 


30 


18 


22 


1 


22 


9 


17 


1 


17 
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Table 4. Unbiased Burst Separation 

AEp< 400 keV AEp> 400 keV 

td 7^0 10 7 

td = 8 12 



Note. — This table shows httle evidence 
that bursts in which emission is delayed evolve 
differently than those which emit promptly. 
However, although the criteria used to divide the 
sample is unbiased, it may not be appropriate 
since most bursts emit very soon after the trigger. 
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Table 5. Bursts Divided According to Emission Time 



AT 


Number 


Min(AEp) 
(keV) 


Med(AEp) 
(keV) 


Max(A^p) 
(keV) 


>0.25 
< 0.25 


10 

27 


80 
120 


300 

525 


950 

2075 


>0.5 
<0.5 


8 
29 


80 
120 


300 

525 


950 
2075 


>0.9 
<0.9 


7 
30 


80 
120 


300 
525 


950 
2075 


> 1 
< 1 


5 
32 


80 
120 


300 

440 


950 
2075 



Note. — This table summarizes the differences in 
evolution of Ep for bursts grouped according to the time 
at which they emit. 
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This preprint was prepared with the AAS I^TJtjX macros v3.0. 

Fig. 1. — Ep as determined from spectra at the finest available time resolution (crosses) and 
from pairs of spectra averaged together (diamonds). The data are from burst 2B910807 and 
the size of the symbols represents the time intervals over which spectra were accumulated 
and the error in Ep. 
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Fig. 2. — Ep as determined from fits to the same spectra with different binnings for burst 
3B930916. The crosses represent values from the unrebinned data while the diamonds 
correspond to the values from spectra which have been rebinned into eight broad energy 
bins. 



Fig. 3. — Ep as measured by detectors at different gains for burst 2B9 10503. The high gain 
detector (crosses) was fit from 25.7-2977 keV and the low gain detector (diamonds) covers 
308-26,214 keV. The high gain detector has better time resolution since it had a stronger 
signal. 



Fig. 4. — Ep as determined from two detectors viewing the burst 2B921207. The high gain 
detector (crosses) covered the energy range from 20.2-1283 keV and the low gain detector 
(diamonds) ranged from 43.1-3300 keV. 



Fig. 5. — The evolution of Ep and the count rate for all bursts in the sample. The histogram 
is the count rate for individual SHERB spectra and the diamonds are Ep. The height of 
the diamonds represent 1-a error bars and their width indicates the accumulation time of 
the SHERB spectra averaged together to achieve the S/N criterion outlined in §3. Note: 
3B931204 extends over two frames. 



Fig. 6. — A comparison of the peak in Ep between intensity pulses within the same burst. 
The diagonal line indicates equality. Harder early pulses lie below this line and harder late 
pulses above it. The bursts used are those for which a +, — , or answer could be determined 
for column 4 in Table 2. Several points line up vertically and horizontally since early pulses 
were individually compared to all those which followed. 



Fig. 7. — The range in Ep for a given emission time for all bursts in the sample. The vertical 
bars give the range in Ep (max-min) and the horizontal bars show the period during which 
significant emission occurred relative to trigger. The crosses mark median values. 



